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ABSTRACT 

We present phase-resolved spectroscopy of two new short period low accre- 
tion rate magnetic binaries, SDSSJ125044.42-I-154957.3 (P„b = 86 min) and 
SDSSJ151415. 65-^074446. 5 (P„b = 89 min). Both systems were previously identified 
as magnetic white dwarfs from the Zeeman splitting of the Balmer absorption lines in 
their optical spectra. Their spectral energy distributions exhibit a large near-infrared 
excess, which we interpret as a combination of cyclotron emission and possibly a late 
type companion star. No absorption features from the companion are seen in our op- 
tical spectra. We derive the orbital periods from a narrow, variable Ha emission line 
which we show to originate on the companion star. The high radial velocity amplitude 
measured in both systems suggests a high orbital inclination, but we find no evidence 
for eclipses in our data. The two new systems resemble the polar EF Eri in its pro- 
longed low state and also SDSSJ121209. 31+013627. 7, a known magnetic white dwarf 
plus possible brown dwarf binary, which was also recovered by our method. 

Key words: stars: binaries: close - stars: magnetic field, white dwarf- stars: low-mass 
- stars: individual: SDSSJ125044.42-hl54957.3, SDSSJ151415. 65-1-074446.5 



1 INTRODUCTION 

Magnetic fields are conunonly observed in white dwarfs. 
Field strengths of 1 MG or greater are found in at least 
10 per cent of isolated white dwarfs ( [Liebert et al.||2003[ ), 
but their formation mechanism is not yet understood. 

It is generally thought that the magnetic fields in white 
dwarfs are 'fossil fields', i.e. a remnant product of the evo- 
lution of chemically peculiar stars of spectral type A and B 
( Moss|1987{|Braithwaite fc Spruit|2004[ ). Afihough a plausi- 
ble formation model, the current space density of Ap and Bp 
stars cannot account for all the magnetic white dwarfs ob- 
served ( Kawka &: Vennes|2004 Wickramasinghe & Ferrario 
|2005[ ). If the progenitors of magnetic white dwarfs are single 
stars, as suggested by this model, we might expect the same 
fraction of magnetic white dwarfs among binary systems as 
found among single stars. Observations, however, show quite 
the opposite. Among the 2248 detached white dwarf plus 
main sequence binary systems found in the Sloan Digital 
Sky Survey (SDSS), not one was found to host a magnetic 
white dwarf ( Rebassa-Mansergas et al. 2011| but see also 



Liebert et al.|2005 where this problem was first addressed) 



Semi-detached magnetic white dwarf binaries, on the 
other hand, are commonly found. They are known as polars 
or intermediate polars, based on the strength of the white 
dwarf magnetic field. They are a subclass of the cataclysmic 
variable stars (CVs), in which mass transfer occurs from an 
M-dwarf companion star to a white dwarf via Roche lobe 
overflow. In polars, the magnetic field of the white dwarf is 
strong enough [B > 10 MG) that it channels the accretion 
fiow directly onto the white dwarf and no accretion disc is 
formed. Magnetic CVs account for ~ 25 per cent of the 
926 known CVs (Ritte r fc Kolb] 12003; v. 7.16, 2011) — a 
considerably higher fraction than is observed among isolated 
white dwarfs. 

This apparent propensity for magnetic white dwarfs to 
form as part of interacting binary systems, has led |Tout et al.| 
( |2008 1 to develop an alternative model for their formation, 
proposing that all high field magnetic white dwarfs have a 
binary origin. The model is based on the standard common 
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envelope (CE) binary evolution theory of Paczynski ( 1976 1. 
When the more massive star in the binary evolves through 
the red giant phase, it expands and envelops both stars. The 
stars continue spiralling around one another inside the CE. 
According to their model, strong dynamo fields are gener- 
ated by the orbiting stars in the core of the CE (see also 
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Potter fc ToutpOlO l. As the stars lose angular momentum 
and spiral closer together into the centre of the CE, the 
magnetic field is frozen into the degenerate stellar core. The 
product of the CE evolution depends on the final separation 
of the components when the envelope is expelled. Binaries 
with a wide separation will have no or very little frozen- 
in magnetic field on the white dwarf, and so leave the CE 
as non-magnetic pre-CVs or non-interacting non-magnetic 
wide binaries. For a strong magnetic field to form in the 
white dwarf, the binary must leave the CE with a small or- 
bital separation. As a result, all magnetic binaries are close 
enough to be interacting, either through Roche lobe over- 
flow, or by wind accretion. In the case where the two stars 
merge inside the CE, it will form a single, massive magnetic 
white dwarf (see also Nordhaus et al.|2011 l. 

Wind-accreting magnetic binaries, as predicted by this 
model, are a recent discovery. The first such system, 



WXLMi (originally known as HS1023-I-3900, Reimers et al 



19991, was discovered in the Hamburg Quasar Survey (HQS, 
Hagen et al.|[l995i) during follow-up spectroscopy of quasar 
candidates. A further eight systems have been discovered 



since (see Schmidt et al. 2007 for a summary). All are 
faint or undetected in X-rays, have cool white dwarfs (TwdJ; 
10 000 K), and display strong cyclotron emission features in 
their optical or near-infrared (NIR) spectra. The accretion 
rates inferred from these characteristics are two to three or- 
ders of magnitude lower than typically observed in polars, so 
initially they were referred to as "low accretion rate polars" 
i Schwope et al.|[2002b |. However, it was later realised that 
the secondary stars in these systems are underfilling their 
Roche lobes (e.g. Schmidt et al.|2005 l, distinguishing them 
from true polars. The inferred low accretion rates are con- 
sistent with magnetic siphoning of the secondary's stellar 
wind (Webbink & Wickramasinghc 20051. These low accre- 
tion rate systems are now recognised as pre-polars, i.e. the 
progenitors of magnetic CVs, and they are expected to start 
Roche lobe overflow in the future as they evolve to shorter 
periods. The nine systems discovered so far lead to an esti- 



mated space density of ~ 10 pc (Webbink & Wickra- 



masinghe|2005 I, more than a factor of ten lower than that of 



polars ( Araujo-Betancor et al.|2005 Thomas & Beuermann 
19981) 



Identifying a newly discovered system as a pre-polar is 
not straightforward, as they share many observable prop- 
erties with low state polars. These "low states" are a well- 
known observed property of polars during which the Roche 
lobe overflow is temporarily interrupted and the accretion 
rate drops by an order of magnitude or more. These low ac- 
cretion states are irregular and unpredictable, and can last 
for days to months ( Kafka fc Hoard|2009 1 . During this time, 
polars appear essentially detached, so it is not immediately 
obvious to which class a newly discovered system belongs. 



For example, EQ Cet was discovered in a low state ( Schwope 
et al.|1999a l and thought to be a pre-polar, but it was later 
detected in a high state ( Schwope et al.|2002a I and dropped 
from the class. 

The reason for the interrupted accretion flow is not yet 
fully understood. Starspots on the secondary star moving 
past the LI point has been suggested as a possible cause 



of this behaviour (Livio & Pringic 1994), but it seems an 
unlikely explanation for the extended low states seen in some 
polars, which can persist for decades. EF Eri is the best- 



known example of such a low state polar. It was known as a 
conventional polar for almost two decades, before it dropped 
into a low state in 1997 ( Wheatley fc Ramsay|1998 |. Apart 
from short bursts of accretion ( Howell et al.|2006 l, it has not 
fully returned to a high state since. With an orbital period 
of 81 minutes, it is close to the orbital period minimum for 



CVs (Gansicke et al._2009|), and a good candidate for having 



a substellar companion ( Beuermann et al.|2000 Schwope & 
Christensen 2010 1. The companion star has not been directly 



detected yet, and a stellar companion cannot be ruled out 
completely ( |HoweU et al.|2006 |. 

Evidence for a substellar companion to a cool magnetic 
white dwarf is also seen in SDSSJ121209. 31-1-013627.7 (here- 
after SDSSJ1212-I-0136). It was first identified as a binary by 
Schmidt et al. (20051, from its optical spectrum. In addition 



to the hallmarks of a magnetic white dwarf, the spectrum 
also showed a narrow Ha emission line, variable both in 
strength and radial velocity on a period of ~ 90 minutes. 
[Schmidt et al.| ( |2005[) favour ed a pre-polar interpretation of 
the system, but |Debes et al. ( 2006| ) showed that the K band 
variability was due to cyclotron emission, suggesting that 
there is ongoing low-level accretion in this system. Variabil- 
ity in the optical is pulse-like and consistent with emission 
from a self-eclipsing heated polar cap, rotating into and out 
of view ev ery 88 minutes (|Koen fc Maxted 2006] [Burleigh] 
et al.|2006 Linnell et al.|2010 1 . This type of optical variabil- 
ity is commonly seen in low state polars, including EF Eri 
(e.g. [Ferrario et al.|1992"{|Howen etal.|2006| ). [Burleigh et al.| 
( 2006 1 also detected X-ray emission from this heated po- 



lar cap, and were able to constrain the accretion rate in 
SDSSJ1212-H0136 to M ~ lO""M0yr-\ Cyclotron emis- 
sion dominates the NIR part of the spectrum and compli- 
cates a direct spectroscopic detection of the companion star, 
but [Farihi et al.[ ( [2008| showed that this part of the spectrum 
can be described by a combination of cyclotron emission and 
an L8 brown dwarf companion star. 

In this paper we present time-resolved optical 
spectroscopy of two more such short period mag- 
netic binary systems, SDSSJ125044.42+154957.3 and 
SDSSJ151415. 65+074446. 5 (hereafter SDSSJ1250-H1549 
and SDSSJ1514-f0744). Both systems have previously 
been identified as magnetic white dwarfs from their SDSS 
spectra ( Vanlandingham et al. 2005 Kiilebi et al. 20091 



and SDSSJ1250-I-1549 was noted by [Steele et al.[ ( [2011[ ) as 
having an infrared excess, possibly due to an unresolved 
companion. The observation that SDSSJ1514-I-0744 is a 
binary system is a new discovery. In the following sections, 
we derive constraints on the binary parameters of these two 
systems and discuss their evolutionary state as polars or 
pre-polars in comparison to similar systems. 



2 DISCOVERY AND FOLLOW-UP 
OBSERVATIONS 

2.1 SDSS and UKIDSS cross-matching 

The two new magnetic binaries were discovered as part of a 
larger project of identifying hydrogen dominated (DA) white 
dwarfs from their colours, rather than through spectroscopy. 
Girven et al. (20111 applied colour cuts to the SDSS Data 



Release 7 (DR7) photometry to select DA white dwarfs, and 



© 2012 RAS, MNRAS 000, ??-[T5l 



Magnetic white dwarf binaries 3 



then fit the ugri magnitudes with a grid of non-magnetic 
white dwarf models to determine the white dwarf parame- 
ters. SDSS spectroscopy is available for 44% of the sources 
which were identified as DA white dwarfs in this way, and 
served as an robust test of the efficiency of the selection 
method. The resulting catalogue of objects was then cross- 
correlated with the Large Area Survey (LAS) of the United 
Kingdom InfraRed Telescope (UKIRT) Infrared Deep Sky 
Survey (UKIDSS) Data Release 8. This allowed us to iden- 
tify white dwarfs with low mass stellar companions, substel- 
lar companions or debris discs, through a near-infrared flux 
excess over the extrapolated best-fit white dwarf model. 

Full details of the cross-matching, selection criteria, 
and the efficiency of the method may be found in |Gir-| 
ven et al. (20111, but we note here that 51 of the 4636 
spectroscopically confirmed white dwarfs were found to 
have an infrared excess. Three of these are of interest 
in this paper: SDSSJ1212-f0136, SDSSJ1250-I-1549 and 
SDSSJ1514-(-0744. These three white dwarfs are clearly mag- 
netic, and were previously identified as such by |Schmidt 



et al. (20031, Vanlandingham et al. (20051 and Kiilebi et al 



(20091, respectively. The SDSS spectra are shown in Fig- 
ure [1 Zeeman splitting of the Balmer absorption lines is 
clearly visible. The dipolar field strengths derived from these 
spectra in their discovery papers are 13, 20 and 36 MG re- 
spectively. 

As discussed in Section [l] SDSSJ1212-I-0136 is a 
known low accretion rate magnetic white dwarf plus brown 
dwarf binary. The two new binaries, SDSSJ1250-fl549 and 
SDSSJ1514-(-0744, are the main subject of this paper. Fig- 
ure |2] shows their spectral energy distributions, as de- 
rived from GALEX [Galaxy Evolution Explorer), SDSS, and 
UKIDSS observations. A model spectrum of a 10 000 K non- 
magnetic DA white dwarf is overplotted on the photometric 
measurements, highlighting the strong infrared excess ob- 
served in both systems. 



2.2 Time-resolved spectroscopy 

SDSS spectra are made up of several subexposures, which 
are individually extracted and combined to produce a final 
spectrum with higher signal to noise. For SDSSJ1250-I-1549 
there are three subspectra available, taken over 56 min- 
utes. For SDSSJ1514-I-0744 there are four subspectra, taken 
over 101 minutes, and another three spectra, of poorer qual- 
ity, taken on a previous day. Inspection of these subspectra 
showed clear radial velocity variations of the Ha emission 
line in SDSSJ1250-I-1549 and a weak indication of a mov- 
ing Ha line in SDSSJ1514-1-0744 (Figure [l]). These results 
suggested that the two new targets may be binaries as well, 
and prompted us to carry out time-resolved spectroscopy on 
them for a more detailed investigation. 

We observed SDSSJ1250-M549 and SDSSJ1514-h0744 
on the nights of 2011 March 24 - 26 using the F0RS2 spec- 
trograph mounted on Unit Telescope 1 of the Very Large 
Telescope (VLT) at Paranal, Chile. Conditions were clear 
with subarcsecond seeing on all three nights. We used the 
GRIS1200R-I-76 grism with a 1" slit and standard 2x2 bin- 
ning, which provided spectral coverage of the wavelength 
range 5872 - 7370A dispersed at 0.73A per binned pixel. 
The resolution was estimated to be 1.7A, based on the mea- 
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Figure 1. Main panel: The SDSS spectra of SDSSJ1212+0136, 
SDSSJ1250+1549 and SDSSJ1514-I-0744, showing Zeeman split- 
ting of the hydrogen absorption lines in dipolar fields of 13, 20 
and 36 MG respectively. All three show weak Ha emission but no 
evidence of a companion star. Small panels: SDSS subspectra of 
the three targets, in the same vertical order. The spectra are cen- 
tred on Ho and are offset for clarity. Variation in radial velocity 
and strength of the Ha line is clearly visible. 
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Figure 2. The Spectral Energy Distribution of SDSS J1250-I-1549 
(top) and SDSSJ1514+0744 (bottom). The measured GALEX, 
SDSS ugriz and UKIDSS YJHK photometry are plotted as red 
dots, the SDSS spectrum as a grey line and a 10 000 K non- 
magnetic DA white dwarf model as a solid black line. The white 
dwarf temperature was derived from a fit to the SDSS ugri pho- 
tometry (see |Girven et aT]|2011| for details). These systems were 
identified as binaries during a search for DA white dwarfs in the 
SDSS, when the cross-match with UKIDSS objects revealed their 
strong infrared excesses. 
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sured full width at half niaxinium (FWHM) of strong sky- 
lines. 

We obtained 24 spectra of SDSSJ1250+1549 on the first 
night, using an exposure time of 360 s, a further six spectra 
on the second night and two on the third night. We used a 
540 s exposure for SDSSJ1514+0744, obtaining six spectra 
on the first night, 13 on the second and 30 on the third. 

The data reduction and spectral extraction were carried 
out using PAMELA ( Marsh|[lj89 1 , which is now available as 
part of the latest STARLlNpfM distribution. The wavelength 
scale was derived from a fifth order polynomial fit to 25 ar- 
clines using MOLL\r| The root mean square (RMS) of the 
fit residuals was 0.012A. The helium-argon-neon arc lamp 
exposures were taken during the day, so to account for small 
shifts due to fiexure of the spectrograph, we shifted the indi- 
vidual spectra by the amount required to move strong sky- 
lines to their known wavelengths. The shifts were typically 
0.1-0.2 pixels. 

The spectra of the two objects display similar character- 
istics: large radial velocity variations, variable Ha emission 
and Zeeman splitting of the white dwarf Ha absorption line. 
No features of a companion star are visible in our spectra. 



2.3 GALEX ultraviolet fluxes 

Both targets are listed in the GALEX all sky survey cata- 
logue, with strong near-ultraviolet {NUV , A^g ~ 2267A) 
and slightly weaker far-ultraviolet (FUV , Agg ~ 1516A) 
fluxes. The ultraviolet magnitudes from GALEX Data Re- 
lease 6 (GR6) are shown in Table [I] and Figure [2] 

2.4 ROSAT non-detection 

SDSSJ1250-I-1549 falls in the field of view of three ROSAT 
{Rontgensatellit) Position Sensitive Proportional Counter 
(PSPC) pointings, the longest of which had an effective ex- 
posure time of 8857 seconds at the position of our target 
(dataset ID RP800393A01, 1992 December 17). It is not 
detected in any of the images. There were also two PSPC 
pointings which included SDSSJ1514-I-0744 in the field of 
view, the longest lasting for 1182 seconds, corrected for vi- 
gnetting effects (dataset ID RP800275N00, 1992 August 12). 
It is also undetected. These non-detections limit the X-ray 
count rates to < 10~^ counts s~^ at the time of the ROSAT 
observations. 



2.5 Optical photometry 

We carried out simple differential photometry (using iraf) 
on our R band acquisition images, and converted the mag- 
nitudes to the SDSS-r band using the transformation equa- 
tions of [Jester et al.|([2005). The magnitudes were constant 
within photometric error over the three nights of observa- 
tions and also consistent with the SDSS DR7 magnitudes as 
listed in Table [l] 

We also show the (unfiltered) lightcurves 
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Figure 3. CRTS lightcurves of siiort period magnetic bi- 
naries, as labelled. SDSSJ1212-(-0136, SDSSJ1250+1549 and 
SDSSJ1514+0744 display no large scale variability, eclipses or 
high states. The lightcurve of EF Eri shows a brief high state in 
2008 November. The short period polar EQ Cet displays large 
amplitude variability on timescales of years. SDSSJ1031-I-2028 is 
currently classified as a pre-polar. Its 42 MG magnetic field shifts 
the cyclotron harmonics to optical wavelengths, so the variability 
in the early part of the lightcurve could be due to variable cy- 
clotron emission. The vertical scale is the same in all the panels 
to directly compare the variability amplitudes. 



of SDSSJ1212-I-0136, SDSSJ1250-H549 and 

SDSSJ1514-I-0744, as obtaine d by the Catalin a Real- 
time Transient Survey (CRTS, [Drake et aL||2009|), in the 



^ http://starlink.jach.hawaii.edu/starlink. PAMELA has been in- 
cluded since the 'Hawaiki' release. 
2 MOLLY was written by TRM and 



available from 



http://www.warwick.ac.uk/go/trmarsh/software/ 



left hand panels of Figure [3| The lightcurves span approx- 
imately three years, with four exposures per visit, taken 
over 30 minutes. The median separation between subse- 
quent visits is 15 days for SDSSJ1250-I-1549 and 18 days for 
SDSSJ1514-I-0744. We find no large scale variability, eclipses 
or high states of these binaries. Figure [3| also includes the 
lightcurves of the low state polar EF Eri, the pre-polar 
EQ Cet and the possible pre-polar SDSSJ1031+2028, for 
comparison. These systems will be discussed in more detail 
in Section [4| 



3 ANALYSIS AND RESULTS 

3.1 Orbital periods 

We measured the radial velocities by fitting a single Gaus- 
sian function to the narrow Ha emission line in each spec- 
trum. The local continuum was accounted for by including 
a pivoting straight line in the fit. We carried out a Lomb- 
Scargle ( Lomb[[1976 Scargle[[1982 1 and orthogonal polyno- 
mial (ORT, Schwarzenberg-Czerny[ [l996[ ) analysis on our 
time series data to search for periodic signals. The ampli- 
tude spectra are shown in Figures [4| and [5[ as these have 
the advantage that the variability at different frequencies 
can be compared directly in physical units. The strongest 
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Table 1. GALEX, SDSS ugriz and UKIDSS YJHK magnitudes of the new magnetic binaries. 



SDSS name 


FUV 


NUV 


u 


9 


r 


i 


z 


Y 


J 


H 


K 


SDSS J125044.42+154957.3 
SDSS J151415. 65+074446. 5 


21.15 
21.64 


19.22 
19.39 


18.65 
19.16 


18.27 
18.84 


18.33 

18.88 


18.36 
18.99 


18.20 

18.88 


17.53 
18.22 


16.65 

18.44 


16.29 
17.33 


14.80 
16.57 



signals are found at 86.3 minutes (SDSSJ1250+1549) and 
88.7 minutes (SDSSJ1514+0744), and are labelled on the 
plots. To confirm that these are the correct periods to as- 
sociate with the orbital period of each system, we compare 
the folded radial velocity curves in Figure [6] A sinusoidal fit 
of the form 



l/(t) = 7 + A'2 sin [27r(i - HJDo)/Po.b] 



(1) 



is overplotted on the measured radial velocities. HJDo is se- 
lected such that it corresponds to the red to blue crossing 
of the radial velocities. In both cases the strongest signal 
results in a radial velocity curve which is well described by 
a sinusoidal fit (Fig. [g] top panel), while folding the radial 
velocities on the one-day cycle count aliases (bottom two 
panels) leads to much larger scatter and x^ values which are 
at least 15 times larger. We also carried out a bootstrap anal- 
ysis to assess the importance of individual measurements in 
determining the orbital period. For each system, we calcu- 
lated the periodogram from 10 000 randomly selected sub- 
sets of the radial velocity curve, and recorded the frequency 
of the strongest peak in each case. For SDSSJ1514-I-0744, 
99.95 per cent of the subsets returned 88.7 minutes as the 
strongest period, and in the case of SDSSJ1250-I-1549, 98.13 
per cent selected the 86.3 minute alias as the strongest sig- 
nal. We are therefore confident that the strongest peaks in 
Figures[4]and|5]correspond to the respective orbital frequen- 
cies. 

To illustrate the phase coverage of each night's observa- 
tions, we plot the folded radial velocities in different colours 
in the middle panels of Figures |4] and |5] The parameters 
of the best-fit circular orbit (Eqn. IT]) and their uncertain- 
ties are shown in Table [2] The deviation from the sine wave 
at phases ~ 0.10 - 0.15 of SDSSJ1514-F0744 is real and 
due to an additional component to the Ha emission (see 
Section 3.3.21. These higher radial velocities are seen at the 



same phase in subsequent orbits and on difTerent nights, and 
were excluded from the sine fit. 

The Zeeman components of the white dwarf's absorp- 
tion lines are clearly variable as well, both in position and 
in strength. By eye, these variations do appear to track the 
orbital motion, but, due to the presence of the rotational 
magnetic field, which modulates the field strength over the 
visible hemisphere, the absorption lines are unreliable as an 
indicator of the white dwarf's radial velocity. 



3.2 White dwarf properties 

SDSSJ1250-H549 and SDSSJ1514-h0744 were both previ- 
ously identified as magnetic white dwarfs from their SDSS 



spectra ( Vanlandingham et al. 2005 Kiilebi et al. 2009 1 



Using a grid of magnetic white dwarf model spectra, along 



with magnetic dipole modelling, Kiilebi et al. ( 2009 1 esti- 
mate their dipole magnetic field strengths to be 20 MG 
and 36 MG respectively. The white dwarf temperature was 
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Figure 4. Top panel: Amplitude spectrum calculated from the 
radial velocities of SDSSJ1250-I-1549. The adopted orbital fre- 
quency is indicated by an arrow. Middle panel: The measured Ha 
radial velocities folded on a period of 86.3 minutes. Data from the 
first night are plotted in black filled dots, the second night in red 
and the third night in green. Error bars are plotted but in most 
cases are smaller than the plot symbols. The dashed line is the 
best-fit sinusoid to the radial velocities (Eqn. [ij. Its parameters 
are listed in Table[2] Bottom panel: Equivalent width variation of 
the emission line, folded on the same ephemeris. 



Table 2. Properties of the new magnetic binaries. 



SDSSJ1250-(-1549 SDSSJ1514-I-0744 Ref. 



Po,b (d) 




0.05995(1) 


0.061610(1) 


1 


HJDo 




2455646.7056(1) 


2455646.89093(3) 


1 


Kern 




360.89 ±2.54 


362.84 ±0.70 


1 


7 (km/s) 
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Figure 5. As in Fig.|4] but for SDSSJ1514+0744. The Ha radial 
velocities (middle panel) and equivalent width (bottom panel) 
are folded on the orbital period of 88.7 minutes. Error bars are 
plotted but in most cases are smaller than the plot symbols. 
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Figure 6. Radial velocities folded on competing one day aliases, 
as labelled. The strongest alias in each case (top panels) is well 
described by a sinusoidal fit, while folding the radial velocities 
on the nearest aliases lead to much larger scatter. Error bars are 
plotted but are smaller than the plot symbols. Points in gray were 
not included in the fit. 



found to be 10 000 K in both systems, using a grid step- 
size of 3000 K and assuming a canonical surface gravity of 
log g = 8.0. The same temperature was obtained by |Girven| 
|et al.| ( |2011[ ) , by fitting non-magnetic white dwarf models to 
the SDSS ugri photometry. The best-fit models are shown 
in Figure [2] We summarise all derived properties of the two 
new magnetic binaries in Table [2] 

For a O.6Af0 white dwarf which has never accreted, a 



tempera ture of 10 00 K corresponds to a cooling age of 
0.6 Gyr (Holberg & Bcrgeroii2006r| A more massive white 



dwarf of the same temperature will be older, increasing to 
1 Gyr for O.SMq , or 2 Gyr for a \Mq white dwarf. The 
surface layers of an accreting white dwarf will be kept hotter 



by compressional heating of the accretion activity ( Townsley 



^ 



Bildsten^^2002J , so the cooling age sets a lower limit on 
the age of the binary. 

We can use the white dwarf temperature and apparent 
magnitudes to estimate the distance to these binaries. As- 
suming a white dwarf radius of O.Oli?0 we calculate the dis- 
tance to SDSSJ1250-f 1549 to be 163±33 pc and the distance 
to SDSSJ1514-I-0744 to be 211 ±42 pc. The uncertainties on 
these distance estimates allow for a ±1000 K variation of the 
white dwarf temperature ( Girven et al.|[20T7 l and a ±0.25 
variation in logg. 

The ROSAT non-detections places an upper limit of 
< 6 X 10"" ergs cm"^ s"^on the 0.1-2.4 keV X-ray fiux, 
which for the assumed distances translate to an accretion 
rate lower than 5 x 1O~^*M0 yr~^. Note however that this 
estimate does not take the accretion luminosity emitted as 
cyclotron emission into account. 



3.3 Location of the Hq emission 

3.3.1 Ha equivalent width 

The high radial velocity amplitude of the Ha emission line, 
as well as the fact that its strength is modulated on the or- 
bital period, suggest that it is emitted from the surface of 
the companion star. To test this hypothesis, we measured 
the equivalent width of the emission line and folded it on 
the orbital ephemeris. To make an accurate measurement, 
we had to account for radial velocity motion of the emission 
line through the absorption trough of the white dwarf. We 
fitted the absorption line with the sum of a narrow and broad 
Gaussian, which was found to describe its shape well. The 
local continuum was normalised using the best fit double 
Gaussian for each spectrum before measuring the equivalent 
width of the emission line in the standard way. The results 
are shown in the bottom panels of Figures H] and [51 folded 
on the orbital period. In both systems the equivalent width 
peaks at phase <j!> = 0.5. As defined by our spectroscopic 
ephemeris this phase corresponds to superior conjunction of 
the companion star, when we are viewing its irradiated side 
face-on. At first glance this type of modulation of the emis- 
sion line strength suggests that the Ha emission is the result 
of irradiation by the white dwarf. However, Ha can also be 
produced by chromospheric activity of the companion star, 
and the two effects may be difficult to disentangle. During 



2007 May and June Howell et al. ( 2008 1 found the Ha emis- 



sion of SDSSJ1212±0136 essentially absent from the optical 
spectrum, over the full orbit. It had returned by 2008 Febru- 
ary, but without a significant change in the overall fiux level 
of the binary or the white dwarf dominated part of the spec- 
trum. The authors therefore conclude that the Ha emission 
cannot be produced by irradiation. Interestingly, the same 



behaviour has been observed in EF Eri as well ( Howell et al 



■^ We used their updated grid of DA white dwarf cooling models, 

available at 

http://www. astro. umontreal.ca/~bergeron/CoolingModels/| 
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20061). It is worth noting that irradiation by a 10 000 K white 



dwarf will only raise the temperature of its companion's ir- 
radiated hemisphere by a few hundred Kelvin, for the orbital 
separations of these binaries, even assuming synchronised ro- 
tation and zero albedo. Using irradiated atmosphere models 
and observations of the polar VV Pup (Porb = 100 min, 
rwd=12 000 K;|Araujo-Betancor et al.|2005| in its low state, 
Mason et al. ( 2008 1 show that irradiation cannot explain the 



SDSSJ1250 + 1544-Ha 
800F 



observed line strengths for white dwarf temperatures below 
20 000 K, and that most of the line emission must have its 
origin in stellar activity. If there is accretion taking place 
in the binary, even at a low level, the accretion spots on 
the white dwarf can easily reach the temperatures required 
to produce ionising photons and cause observable irradia- 
tion of the companion star, even if the white dwarf itself 
is too cool. The shape of the phase-folded equivalent width 
lightcurves of SDSSJ1250+1549 and SDSSJ1514+0744 does 
mean that if the Ha emission here is due to stellar activity, 
it must be restricted to, or greatly enhanced, on the hemi- 
sphere of the companion facing the white dwarf only. Such 
a scenario could result from magnetic interaction between 
the white dwarf and its companion, which would induce ac- 
tivity in the companion preferentially on the white dwarf 
facing hemisphere ( [Howell et al. 12006] ). 

3.3.2 Doppler tomography 

Doppler tomography ( [Marsh fc Horne||1988l [Marsh[[2005l ) 
provides an alternative way to locate the source of the emis- 
sion. The method uses the fact that the observed line profile 
at a given orbital phase is the result of the integrated emis- 
sion at the line wavelength, moving at different velocities 
in different parts of the binary system. The binary rotation 
allows us to view the system from different angles and sepa- 
rate the different velocity components. Such a reconstructed 
"velocity image" can then be compared to the expected ve- 
locities of the binary components to infer the origin of the 
emission. 

Our Doppler tomograms are shown in Figure (tI along 
with the trailed spectra they were computed from. Here we 
simply subtracted polynomial fits to the continuum from 
the spectra without removing the white dwarf absorption or 
any other features of the spectrum. For comparison, Doppler 
maps are overplotted with the Roche lobe of a secondary 
star in a binary with mass ratio q = M2/M1 = 0.2, for an 
assumed minimum and maximum value of the A'-correction 



(see Section 3.3.31. The corresponding velocity projections 



of a ballistic stream emanating from the inner Lagrangian 
point are also shown. 

The strong emission line dominating the trailed spec- 
tra maps to a single bright point in the tomogram, con- 
sistent with the Ha emission originating on the compan- 
ion star. No other contribution is seen in the tomogram of 
SDSSJ1250-hl549 (left hand panel), but SDSSJ1514-h0744 
displays additionally some faint emission which appears to 
be active mass transfer between the two stars. The stream 
does not follow a ballistic path for any realistic set of binary 
parameters. It has a lower velocity than the ballistic stream 
in all cases. This feature could be evidence for wind accre- 



tion by the white dwarf ( Webbink & Wickramasinghe 2005 1. 
Such an interpretation of the tomogram should be treated 
with some caution. A magnetically controlled outflow will 
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Figure 7. Doppler tomography and trailed spectra of 
SDSSJ1250+1549 (left) and SDSSJ1514+0744 (right). The Roche 
lobe of the secondary star in a binary with q = 0.2 is overplotted 
on the maps as a reference, for two values of the iC-correction. 
The smaller Roche lobe (centre of mass marked by a X symbol) 
has a zero A'-correction, i.e. assumes that the Ha emission is dis- 
tributed uniformly over the surface of the star. The larger Roche 
lobe (star symbol) is calculated assuming the emission arise from 
the irradiated inner face of the secondary star only, with no con- 
tribution from the cooler opposite hemisphere. The + symbol at 
(0,0) indicates the centre of mass of the system. The predicted 
trajectory of a ballistic gas stream emerging from the inner La- 
grangian point in each case is shown as well. The trailed spectra 
in the lower panels are phase binned, and two cycles are plotted 
for clarity. Note the faint S-wave of the mass transfer stream in 
the trail of SDSSJ1514-I-0744. 



have a velocity component out of the plane of the binary 
orbit, violating a key assumption of Doppler tomography. 
Nevertheless, this technique has been successfully applied 
to polars with magnetically channelled accretion flows in 
the past (e.g. Schwope et al. 1999b I, often showing a bal- 



listic gas stream leaving the LI point before it is disrupted 
and channelled by the magnetic field within the Roche lobe 
of the white dwarf. The additional emission component is 
also visible as a faint S-wave in the trailed spectra and as a 
radial velocity disturbance in Figure [5] 



3.3.3 K- correction of the Ha line 



We have shown in Section |3.3| that the Ha emission origi- 
nates mainly from the inner hemisphere of the companion 
star. Therefore, the measured emission line radial velocity 
amplitude (K^^) corresponds to the velocity of the centre of 
light, rather than the centre of mass of the secondary star. 
As the centre of light is closer to the centre of mass of the 
binary, K^ia underestimates the true centre of mass radial 
velocity, K2. 

In order to place a dynamical constraint on the white 
dwarf mass, a correction factor needs to be applied. This cor- 
rection factor is a function of the mass ratio q = AI2/M1 — 
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K1/K2 and the distribution of the Ha emission across the 
surface of the star. The latter is difficult to determine ex- 
actly, but it can be constrained through geometrical consid- 



erations ( Wade & Home 1988 Munoz-Darias et al. 2005 1 



The corrected radial velocity amplitude of the secondary 
star is given by 



K2 = 



A"c. 



1 - (1 + g)/i?2/a 



(2) 



where fR2 is the distance of the centre of light from the 
centre of mass of the secondary star, with 0^/Sjl./ = 
represent the case where Ha is emitted uniformly across the 
surface of the secondary star, so that K2 = -ft^cm, while / = 1 
is the opposite extreme where all the emission originates in a 
small spot on the secondary star closest to the white dwarf. 
For the case where the hemisphere facing the white dwarf 
emits uniformly, and there is no contribution from the cool 
hemisphere, / = 4/37r ^ 0.42 (|Wade fc Horne||1988[). We 



will assume this value as representative in the calculations 
that follow. For two stars separated by a distance a, orbiting 
their common centre of mass with a period Poi b , the orbital 
velocity of the secondary star is given by 

27ra sin i 



K2 



(3) 



where i is the binary inclination. Combining equations [2] and 
[3] allows us to write the measured emission line velocity K^^a 
as a function of the mass ratio q and the binary inclination 



K, 



2-Ka sin i 



l-(l + g: 



fR2 



We return to this equation in Section [4?^ 



(4) 



3.4 Spectral type of the secondary star 

No features of the secondary star, apart from the Ha emis- 
sion, are detected in our time series spectra. The UV-optical 
part of the SED is well described by a 10 000 K white dwarf 
in both systems but there is a large excess in the UKIDSS 
bands (Figure [2|. [Steele et al. (2011 1 modelled the excess in 
SDSSJ1250+1549 with an M8 companion, but were only 
able to match the JH flux with this spectrum. The en- 
hanced K band flux is likely due to cyclotron emission, but 
its overall contribution to the spectrum cannot be deter- 
mined from these data alone. In general, the cyclotron emis- 
sion is strongly modulated on the orbital period, making the 
UKIDSS 'snapshot' observations unsuitable for SED mod- 
elling. 

For SDSSJ1514-I-0744, the cyclotron harmonics wiU ap- 
pear at shorter wavelengths than in SDSSJ1250+1549, due 
to its stronger magnetic field. Cyclotron emission will almost 
certainly contribute to some or all of the NIR flux measure- 
ments but by an (as yet) unknown amount. We note however 
the very low J band flux of SDSSJ1514-h0744 (see Figure[2|. 
There is no error flag on this measurement in UKIDSS, so it 
appears to be real. If so, it offers a constraint on the spectral 
type of the companion. Assuming that the cyclotron humps 
fall outside of the J band in this system, so that the emission 
in this band is from the companion star only, the measured 
flux is consistent with an L3 type star or later at the de- 
rived distance of 211 pc. NIR spectroscopy will be essential 



to model the cyclotron emission and verify this estimate, as 
well as to search for other spectral features of the companion 
stars of both systems. 



4 DISCUSSION 

4.1 Distinguishing between low state polars and 
pre-polars 

The component stars in polars and pre-polars are thought to 
be the same, but in the latter, accretion occurs via magnetic 



siphoning of the secondary's stellar wind ( Webbink & Wick- 
ramasinghe|2005 1 rather than through Roche lobe overflow. 
The inferred low accret ion rates of the pre- polars (M ~ 
10-14 _ 10""Mo yr"MSchmidt et al.||2007| are compati- 



ble with this model, although an accretion stream has been 
detected in at leas t one system, SDSSJ204827.9+005008.9 
( Kafka et al.|2010 |. 

There are striking similarities between the two classes, 
such as the cool white dwarf, low inferred accretion rate and 
magnetic field strengths such that the cyclotron harmon- 
ics are shifted into the optical/NIR part of the spectrum. 
We compare graphically the properties of SDSS J1212-I-0136, 
SDSSJ1250+1549 and SDSSJ1514+0744, with those of 
known polars and pre-polars in Figure [8] We show only the 
polars for which a reliable temperature estimate is available, 
taken from Table 1 of Townslcy & Gansicko ( 2009^ . As was 
already pointed out by |Schwope et al.| ( |2009[ ), the pre-polars 
are generally found at a slightly lower temperature than even 
the coolest polars. This can be understood in terms of ac- 
cretion heating — a higher accretion rate in polars leads to 



more compressional heating of the white dwarf (Townsley 



fc Bildsten 2002 Townsley fc Gansicke 20091, raising the 



observed effective temperature. 

The orbital period is a potential discriminator between 
the two classes. The pre-polars all have periods longer than 
2.5 hours, with the exception of SDSSJ103100. 6+202832. 2 
(shorthand notation: SDSSJ1031-I-2028) at 82 minutes 



(Schmidt et al. 2007 Linnell et al. 20101. It seems plausi- 



ble that the cooler, long period systems are still evolving 
towards Roche lobe contact, but it is important to con- 
tinue monitoring these systems to reveal any high states 
which may prove them to be low state polars instead, as 
happened in the case of EQ Cet ( [Schwope et al.||1999a[ 
,2002a). Apart from being the only short-period pre-polar, 
SDSSJ1031-I-2028 is also the only one without evidence for 
a late type secondary star in its optical spectrum. With a 
magnetic field of 42 MG, the optical emission is dominated 
by cyclotron harmonics, and the spectrum has no clear fea- 
tures which can be attributed to the companion star. In fact, 
the m = 4 harmonic is centred very close to Ha, complicat- 
ing an analysis of the emission line variability. Based on the 
absence of a molecular band at 7600 A (which falls between 
the cyclotron harmonics) Schmidt et al.| ( |2007[ ) are able to 
limit the spectral type of the secondary star in this binary 
to ^M6, meaning that it is underfilling its Roche lobe for 
the 82 minute orbital period. We also note that its magnetic 
field is lower than the other pre-polars, but higher than that 
of EF Eri, SDSSJ1212-I-0136 and the two new binaries dis- 
cussed here. 

It is interesting to note the position of EF Eri in the 
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Teff-Porb (second) panel of Figure [s] It is the shortest pe- 
riod, lowest temperature polar on the plot, very close to 
SDSSJ1031-I-2028. EF Eri was known as a typical high stat e 
polar from the time of its discovery (Williams et al.||1979l, 



until it dropped into a low state in 1997 ( Wheatley & Ram- 
[sayli ggs"). It has remained in a low state ever since, showing 
only brief bright states several years apart, but no return to 
its former high stataj One of these short bright states was 
detected by the CRTS in 2008 November and is shown in 
Figure [3] The lightcurve of EF Eri shows no long timescale 
variability during its low state like the polar EQ Get in 
the panel below it. In fact, it displays greater similarity to 
the lightcurves of SDSSJ1212+0136, SDSSJ1250+1549 and 
SDSSJ1514+0744. The lightcurve of SDSSJ1031+2028 dis- 
plays more variability, likely caused by the variable cyclotron 
emission dominating its optical spectrum. 

The general similarity between EF Eri and 
SDSS J 103 1-1-2028 was also noted by [Schmidt et al.| 



1 20071 and it remains a possibility that SDSSJ1031-I-2028 
could turn out to be a polar in an extended low state. 



4.2 Evolutionary status 

Where do the new binaries fit into this picture? The pa- 
rameters of the three magnetic binaries, SDSSJ1212-I-0136, 
SDSSJ1250-H549 and SDSSJ1514-h0744, are plotted in Fig- 
ure ID in filled black symbols. With their 10 000 K white 
dwarfs, strong magnetic fields and orbital periods 80 — 90 
minutes, they occupy the same parameter space as EF Eri 
and SDSSJ1031-I-2028. Apart from the variable Ha line, no 
spectral features of the companion star have so far been de- 
tected in their spectra. Flux considerations limit the spectral 
type of the companion star in SDSSJ1212+0136 to L5 - L8 
I Schmidt et al.||2005| [Farihi et alT]|2008[ ), but it is not clear 
yet whether this system is a detached binary or whether 
it is a polar currently in an extended low state. Our cur- 
rent data on SDSSJ1250+1549 and SDSSJ1514-h0744 can- 
not constrain the parameters of these systems very tightly, 
so we consider below three possible evolutionary scenarios: 
i) a polar in an extended low state, in which we assume that 
the secondary star fills its Roche lobe, ii) a detached binary 
which has not yet evolved to Roche lobe contact, and for 
which we assume a mass-radius relation for the companion 
star appropriate for isolated late type stars, and in) a polar 
which has evolved past the orbital period minimum of CVs, 
and is now cooling and fading. 



4-2.1 Low state polar 

Under the assumption that these magnetic binaries are po- 
lars, we can take advantage of the fact that the secondary 
stars must be filling their Roche lobes, and approximate 7?2 
by the volume-equivalent Roche lobe radius, Rl- We use the 
expression derived by |Sirotkin fc Kim| ( |2009"l ), which models 
the secondary star as an n = 3/2 polytrope. 



A long term monitoring lightcurve from the SMARTS 
telescope was presented by F. Walter in 2009 March at 
the 14th North American Workshop on Cataclysmic Vari- 
ables and Related Objects. The presentation available at 
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Figure 8. White dwarf magnetic field (top panel), white dwarf 
temperature (second panel) and orbital period distribution (bot- 
tom three panels) of polars, pre-polars and the three magnetic 
binaries discussed in this paper. The light gray shading indi- 
cates the CV period gap. The 'new' magnetic binaries (filled 
black squares, black histogram) have roughly the same period 
and temperature and similar magnetic fields. The pre-polars (gray 
dots, gray histogram) appear preferably at longer orbital peri- 
ods, slightly cooler white dwarf temperatures and higher mag- 
netic fields. The polars (open triangles, white histogram) span a 
wider range in temperature and magnetic field than the other two 
classes. The parameter space Twd^ 10 000 K, Porb ^ 90 min, con- 
tains SDSSJ1212+0136, SDSSJ1250-I-1549, SDSSJ1514+0744, a 
known polar in a low state (EF Eri), another polar previously 
observed in both high and low states (FL Get) and a presumed 
pre-polar (SDSSJ1031-I-2028). 



Rl 



0.5126 g° 



0.6710q07349 4.in(l + q0.: 



(5) 



This approximation to the volume-equivalent radius offers 



an improvement over the Eggleton ( 1983 \ relation for donor 



http://www.noao.edu/meetings/wildstars2/talks/wednesday/walter 



stars below the CV period gap which are fully convective 
and therefore less centrally concentrated. 

Using equation |4] we can calculate all possible {q,i) 
combinations which give the measured emission line veloc- 
ity Ti'cm. The results are shown in Figure [9] We fixed Mi at 
O.6Af0 , O.75M0 and 0.9Mq in three separate calculations, 
and connect the possible (M2,i) pairs for each case with a 
black solid line. The binary inclination can be constrained 
based on the fact that we find no evidence for eclipses in 
our spectra or in the CRTS lightcurves. This restricts the 
inclination to z < 72°. The inclination constraint of course 
depends on the on the mass ratio and on the radius of the 
ppfe3Condary star. The limit we plot in Figure [9] corresponds to 
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Figure 9. The emission line velocity amplitude i^om is a function 
of both binary inclination i and the mass ratio q (see text for 
details). The solid black lines are the K^^ contours corresponding 
to the velocity amplitude measured from our spectra, for each 
pair of M2,i values, and the white dwarf mass M\ as labelled. 
The corresponding error on A'em is indicated as gray shading. 
The area above the dotted horizontal line is ruled out as we find 
no evidence for an eclipse in our data. The binary parameters, 
as derived from the [Knigge et al.| l |2011[ | donor star sequence, is 
indicated by a star symbol. The 'RLF' is a reference to the fact 
that this calculation assumes Roche lobe filling. 



the lowest inclination for the range of masses considered for 
which no eclipse of the secondary star would be observed. 
For the lowest values of M2 eclipses would not be observed 
for inclinations up to 77°. The lower limit on the inclination 
is set by the Chandrasekhar mass limit for the white dwarf, 
restricting it to i > 38° in SDSSJ1250+1549 and i > 39° in 

SDSSJ1514+0744. 

The CV donor star sequence of [Knigge et al.| ( |2011[ ) is 
a semi-empirical sequence describing the physical and pho- 
tometric properties of the secondary stars in cataclysmic 
variables, as a function of orbital period. It was derived us- 
ing non-magnetic CVs only, so strictly it only applies in the 
case of non-magnetic systems, but we include it as it may 
provide some insight into the properties of these binaries. 
The effect of the white dwarf magnetic field is to reduce the 
efficiency of magnetic braking of the secondary star by par- 



tially or fully controlling the mass flow ( Wickramasinghe fc 
Wu|1994 l. The reduced angular momentum loss results in a 
reduced mass transfer rate and therefore donor stars which 
are less out of thermal equilibrium than their non-magnetic 
counterparts at the same Porb- 



For the measured orbital periods of 86 and 89 min- 
utes for SDSSJ1250+1549 and SDSSJ1514-f0744 respec- 
tively, the typical non-magnetic CV donor has a mass of 
0.077Mgand O.O82M0 , with radii O.124i?0and 0.1297?© . 
The calculation assumes a representative white dwarf mass 
of Ml = ^.ToMq , so we indicate these solutions as sin- 
gle points in Figure [9] The corresponding binary inclination 
values are 60° and 62°. 

As a result of the lower mass transfer rate, the white 
dwarfs in polars are generally found to be cooler than 
non-magnetic CV white dwarfs (Sion 2004). A 10 000 K 
white dwarf temperature implies an average accretion rate 
over the thermal timescale of the envelope (~ 10^ year) 
of M = 1.2 X lO""M0yr"^ for a O.OMq white d warf, or 
M = 4.8 X 1O""M0 yr-i for a O.6M0 white dwarf (JTowns- 



ley & Gansicke 2009). These estimates are orders of mag- 



nitude higher than what we infer from the ROSAT non- 
detections, but not atypical of the secular accretion rate of 
polars ( Warner|2003 p. 347). The implication is that, unless 
they are younger than 0.6 Gyr, the white dwarf in these 
binaries must have accreted at a much higher rate in the 
relatively recent past. 

4-2.2 Detached binaries 

In this section we consider these binaries under the assump- 
tion that they are detached pre-polars. In this case, the sec- 
ondary stars have never filled their Roche lobes and the only 
accretion that has taken place would be through magnetic 
capture of the secondary's stellar wind. Wind accr etion pro- 
ceeds at a much lower rate (~ 6 x 10^ ^^Af, 



% yr 



Webbink 



& Wickramasinghc'2005i) than typically seen in systems with 



Roche lobe overflow, so causes very little accretion-induced 
heating of the white dwarf. Unlike for CVs, the white dwarf 
temperature can therefore be used as an approximate age 
indicator. For a O.6M0 white dwarf (resp. O.9M0 ) a tem- 
perature of 10 00 K corresponds to a. cooli ng age of 0.6 Gyr 
(resp. 1.3 Gyr) |Holberg fc Bergeron||2006p . 

Due to this low mass loss rate, the structure of the com- 
panion stars is probably closer to that of non-interacting low 
mass stars, than to that of the companion stars of CVs. In 
Figure [lO] we show the range in inclination and compan- 
ion star mass which can result in the observed emission line 
velocity K^^n- Rather than assuming a Roche lobe filling sec- 
ondary star as in the previous section, we use here the evo- 
lutionary models of Ba raffe et al.||T998[ ), for low mass stars 
of ages 0.6 Gyr and 5 Gyr to determine the mass and radius 
of the companion stars. For simplicity, we show only the re- 
sults for 0.6 Gyr in Figure [lO] The 5 Gyr model results in a 
very similar diagram and does not affect the range of possi- 
ble companion star masses derived below. The mass-radius 
relation for both models are shown in Figure [TT| along with 
measurements of non-interacting low-mass stars, as com- 



piled from literature by Knigge et al. (20111. Although the 
two models are very similar over a large part of the mass 
range under interest, they diverge slightly for masses below 
AI2 < 0.15, with the younger stars having a larger radius for 
a given mass. This slight difference in radius translates to 
a large difference in the mean stellar density for low mass 
stars. The mean density is proportional to M/R^, so for 
a given (low) mass, the older stars are smaller, and there- 
fore much denser. This difference is illustrated in Figure [12] 
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For comparison, we also show the average density of the 



Knigge et al. (20111 CV donor star sequence, as well as a 



360.89 [SDSSJ1250+1549] 



lower limit on the densities of the binaries considered in this 
paper. These limits, as indicated, were calculated from the 
period-density relation, which assumes that the stars are 
Roche lobe filling. If they are instead detached, as we as- 
sume in this section, their radii will be smaller for the same 
period and the densities higher. For non-detached binaries, 
the period-density relation therefore provides a lower limit 
on the mean stellar density. 

Since the orbital periods are well determined, we can 
also use the period-density relation to set an upper limit 
on the mass of the companion star. In Figure |10[ the thick 
vertical dashed line indicates the mass of the companion star 
which have the same density as a Roche lobe filling star at 
the orbital period measured for these binaries. As before, we 



have used the mass-radius relation from Baraffe et al. ( 1998 1. 



It is clear from Figure[T2|that, for the same mass, the density 
of non-interacting stars (dashed and dotted lines) is higher 
than the CV donor stars (solid line). This is a result of the 
mass transfer which drives the donor stars of CVs out of 
thermal equilibrium, causing them to expand slightly. From 



eclipse measurements, Littlefair et al. ( 2008 1 find the donor 
stars in CVs to be ~ 10% larger than isolated M-dwarfs of 
equal mass. 

The mass limit derived this way is a strict upper limit. 
If the companion stars were more massive, their radii would 



be larger (Figure 111 and hence their densities lower (Fig- 
ure 12 1. A larger radius means that the star would exceed the 



size of its Roche lobe so the binary would no longer be de- 
tached. We can constrain the mass further by requiring that, 
for each mass-radius pair, the binary remains detached, i.e. 
that the companion star radius is smaller than the Roche 
lobe radius. This is only true for 0.04 < M2/M0 < 0.12 (in 



both systems), shown in light gray shading in Figure 10 



The mass range derived above allows for both a stellar 
or a substellar detached companion to the white dwarf. But 
how common are such binaries? White dwarf plus M-dwarf 
binaries are plentiful, and a common product of the com- 
mon envelope evolution phase of binary stars. The SDSS 
DR7 catalogue of white dwarf-main sequence binaries, for 



example, contains 2248 entries ( Rebassa-Mansergas et al 



2011), of which approximately one third are post common 



envelope binaries (Schreiber et al. 20101. However none of 



these binaries have a magnetic white dwarf like the binaries 
considered here. Brown dwarf companions to white dwarfs 
are also a rare find. From 4636 spectroscopically confirmed 



white dwarfs in the SDSS DR7, Girven et al. (2011 1 find less 



than 2 per cent to have an infrared excess consistent with 
a brown dwarf companion. In an earlier survey, [Farihi et al.| 
( 2005[ ) find fewer than 0.5 per cent of white dwarfs to have 
brown dwarf companions. Only seven white dwarf-brown 
dwarf binary systems have been spectroscopically confirmed, 



and only two of these are close binaries, GD 1400 ( Farihi & 
|Chri stopher'l2004l) and WD 0137-349 ( [Maxted et al.|2006 r 



Although we cannot rule out that SDSSJ1212-I-0136, 
SDSSJ1250-hl549 and SDSSJ1514-h0744 are detached sys- 
tems, it seems rather unlikely that we would serendipitously 
find three near-identical detached magnetic binaries, all at 
a similar distance, in only ~ 2700 degrees of sky (see Sec- 
tion 4.3 1. The implied space density of such a find is much 




K,,„ = 362.84 [SDSSJ1514+0744] 
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Figure 10. As for Fig.|9] but here we use the mass-radius relation 
of an isolated M-dwarf of age 0.6 Gyr, rather than assuming a 
Roche lobe filling companion star. The thick dashed vertical line 
is a strict upper limit on the mass of the companion, from the 
density of the companion star. Requiring that the system remains 
detached, sets further limits on the mass (see text), leaving only 
the light gray shaded area as possible parameters. 




Figure 11. The mass-radius relation of non-interacting low- 
mass stars. The points with errorbars are measurements from 
literature, compiled by [Knigge et alT]l|2011[|. T he dashed line is 
the 0.6 Gyr isochrone from |Baraffe et al.| ( |1998^ , which is a lower 
limit on the age of the binaries considered here. The dotted line 
is for 5 Gyr. Note the deviation of the two models at low masses. 



higher than suggested by the small number of known pre- 
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Figure 12. Average stellar density for different models of the companion star. The dashed and dotted lines are the 0.6 and 5 Gyr low 
mass stellar models from.Baraffe ct al. (1998) and the solid black line the CV donor sequence from [Knigge et al.| ( |20lH . Overplotted 
as filled black circles are the average densities of CV donor stars with a reliable mass and radius measurement from eclipse modelling 
( [Savoury et al.|2011| l. The open circles are the companion stars in detached white dwarf binaries, the triangles represent low mass stars 
in other types of binaries and the squares are the measurements of single stars. The data were taken from the compilation by Kni gge| 
|et al.| ( |201l'| |. A lower limit on the densities of SDSSJ1031+2028, SDSSJ1212+0136, SDSSJ1250+1549 and SDSSJ1514+0744, calculated 
from the CV period-density relation, are indicated by thin gray lines. The density of EF Eri is also shown, and as a known polar, its 
density must be close to the value shown. 



polars or non-magnetic white dwarf plus brown dwarf bina- 
ries. 

The constraints derived from observable parameters 
leave more scope for these binaries to be Roche lobe fill- 
ing (or nearly so) than detached. We note however that if 
these three binaries are pre-polars, it will extend (to lower 
values) the range of magnetic field strength and orbital pe- 
riod pre-polars are observed to have. 



4-2.3 Post-bounce polar 

In the standard model of CV evolution a system evolves 
to shorter periods through angular momentum loss, until 
a minimum period Pmin, when the companion star can no 
longer sustain hydrogen burning in the core and becomes 
degenerate. The current observational estimate of this min- 
imum period is Pmi„ = 82.4 minutes ( Gansicke et al.|[2009 |. 
Further mass loss from the degenerate companion causes 
the binary orbit to expand and the system to evolve back to 
longer periods. It is observationally challenging to directly 
detect such a degenerate secondary in a CV, especially since 
the optical luminosity is usually dominated by an accretion 
disc. The absence of an accretion disc in polars greatly im- 
proves chances of a direct detection, but even so it remains 
an elusive result. The third evolutionary state we consider 



for these binaries is that they are polars which have evolved 
past the period minimum. 

If we assume that the CV evolutionary track of 



Knigge et al. (20111 applies to SDSSJ1250+1549 and 



SDSSJ1514-I-0744, their orbital periods imply a binary age 
of > 3.3 Gyr. In the case of polars, the efficiency of magnetic 
braking is reduced by the magnetic interaction between the 
two component stars, so the evolution probably proceeds 
more slowly, making the binaries even older than this esti- 
mate. Nevertheless, even at 3 Gyr, the spectral type of the 
donor is expected to be as late as T, with a mass of only 
hh ~ O.O4M0 . The mass transfer rate from the donor drops 
as the system evolves past Pmin, so the low current accretion 
rate may be expected. 

It is interesting to note that for the non-magnetic CVs of 
the Knig ge et al.| ( |2011[ ) evolutionary track, the white dwarf 
temperature of a post-bounce CV at these orbital periods, 
is expected to be in the range 9160-11 170 K. Although the 
observed temperatures of our magnetic binaries fall precisely 
in this range, it is not clear that this estimate will apply to 
magnetic CVs. The reduced angular momentum loss rate in 
magnetic CVs results in a lower accretion rate and less com- 
pressional heating of the white dwarf, so T„d is expected to 
be lower. By the same token, the evolution of magnetic CVs 
proceeds more slowly, so a magnetic CV is likely to be older 
and hence cooler than a non-magnetic CV at the same or- 
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bital period. An isolated O.6Af0 white dwarf will have cooled 
to 5500 K in 3 Gyr (7000 K for a 0.9Mq white dwarf), so 
the observed 10 000 K must be result of relatively recent 
accretion at a high rate. 



4.3 Space density 

The overlap in the SPSS DR7 and UKIDSS /LAS DR8 foot- 
prints is ~ 2700 degrees ( Girven et al. 2011 ). Our crossmatch 



revealed three magnetic white dwarfs with a low mass com- 
panion within this area, suggesting that the space density of 
such systems may be quite high. We estimate the effective 
survey volume searched by [Girven et al.| ( |2011[ ) by integrat- 
ing the volume enclosed by a spherical cone, taking into 
account the exponential drop-off of the space density per- 
pendicular to the galactic disc (for more details, see Sect. 
5.1 of Gansicke et al.|[2009 |. We weigh the contributions of 
individual galactic latitudes by the fractional area covered 
by UKIDSS/LAS DR8. Adopting d = 200 pc as a typical 
distance to the magnetic binaries discussed in this paper, 
and a scale height of either H^ = 190 pc ( Patterson]|1984 1 
or Hz = 260 pc ( [Pretorius et al.||2007[ ), we find an effective 
survey volume of 5.9 x lO'^pc^ and 6.6 x lO^pc'^. Because 
of their intrinsic faintness, all three systems are within one 
scale height, i.e. the uncertainty in H^ does not significantly 
affect our estimates, and we find p ~ 5 x 10~®pc~^. 

This value is substantially higher than the estimated 
space density of "genuine" polars, ~ 1.3 — 1.6 x 10~''pc~^ 



( Thomas &: Beuermann|1998 Araujo-Betancor et al.|2005 l, 
and, in fact, comparable to the space density of non- 



magnetic CVs, 4^^pc ^ as est imated by 



(20071 and Pretorius & Knigge (20121 from an analysis ot 



Pretorius et al. 



two ROSAT X-ray surveys. 

An alternative approach of assessing the relative fre- 
quency of the low-accretion rate magnetic binaries found by 
the SDSS, is to compare their number to that of equally 
bright non-magnetic SDSS CVs whose optical spectrum 
is dominated by their accreting white dw arfs. Within the 
Legacy footprint of SDSS DR7 (8032 deg^,|Abazajian et aL 



2009), Szkody et al. (20111 have identified 15 CVs with 



g ^ 18.84 that clearly exhibit broad Balmer absorption 
lines from the white dwarf photosphere. This gives a sur- 
face density of 1.9 x 10~^deg~^, which is comparable to the 
1.1 X 10~^deg~'^ calculated for the low-accretion rate mag- 
netic binaries. The spectroscopic completeness of SDSS is 



similar for both types of objects (^Girven et al. 2011 1, and be- 



cause of their clear observational hallmarks (non-magnetic 
CVs: Balmer emission lines from the accretion disc; low- 
accretion rate magnetic binaries: infrared excess due to cy- 
clotron emission), they have been identified from their SDSS 
spectra and other obeservations to a high degree of com- 
pleteness as well. 

Unless the sky distribution of the three low-accretion 
rate magnetic binaries discussed in this paper is subject to 
small-number anomalies, we have to conclude that they rep- 
resent an abundant sub-population of close white dwarf bi- 
naries whose existence is not readily explained within the 
current theoretical framework. 



5 SUMMARY AND FUTURE OBSERVATIONS 

We have presented VLT/F0RS2 spectroscopy of two mag- 
netic white dwarfs, SDSSJ1250+1549 and SDSSJ15144-0744 
and have shown that they are binary systems, hosting low 
mass companion stars. We measured their orbital periods 
from a variable Ha emission line in their spectra and showed 
that in both systems this emission originates from the sur- 
face of the companion star. Their orbital periods are 86.3 
and 88.7 minutes respectively. No other features that can 
be attributed to the companion star are seen in our spectra. 
SDSSJ1514-I-0744 additionally shows evidence for a high ve- 
locity Ha component. Its velocity is lower than that of a bal- 
listic stream, and may be interpreted as evidence of magnetic 
siphoning of the companion's stellar wind. This strongly sug- 
gests that there is ongoing accretion in this system. In both 
binaries, the Ha equivalent width is also modulated on the 
orbital period, and peaks when the companion star is at 
superior conjunction. The Ha emission must therefore orig- 
inate mainly from the hemisphere of the companion facing 
the white dwarf, although it is not yet clear whether it is 
the result of irradiation by a hot spot on the white dwarf or 
magnetic activity on the companion star. The large radial 
velocity amplitude of the Ha line (> 360 kms~^) points to 
a high inclination in both systems, but we find no evidence 
of eclipses in the CRTS optical lightcurves of these binaries. 

Our current data cannot constrain the parameters of 
the binary components very tightly, so we have considered 
in detail three possible evolutionary scenarios: polar in an 
extended low state, detached pre-polar, or polar that has 
evolved past the orbital period minimum of cataclysmic vari- 
ables. Although none of these possibilities can be ruled out, 
we favour the low state polar interpretation for the observed 
similarities between these binaries and EF Eri, a well known 
polar which have been in a low state since 1997. 

The properties of these two binaries are also near- 
identical to SDSSJ1212+0136, a known magnetic white 
dwarf plus possible brown dwarf binary. Incidentally, 
SDSSJ1212-I-0136 was also recovered by the method used to 
find these two new binaries. Just like SDSSJ1212+0136, the 
SEDs of SDSSJ1250+1549 and SDSSJ1514-h0744 display a 
large infrared excess, which we have interpreted as cyclotron 
emission from the magnetic field of the white dwarf and pos- 
sibly emission from a late type companion star. 

Clearly these binaries will benefit greatly from further 
observations. In particular, an X-ray detection will give an 
improved estimate of the accretion rate, and time-resolved 
infrared spectroscopy will allow us to model the cyclotron 
emission, which will also constrain the accretion rate. In- 
frared spectroscopy will also allow us to look for spectro- 
scopic signatures of the companion star, or place limits on 
its spectral type from flux considerations. Finally, an im- 
proved distance estimate will allow us to better constrain 
the white dwarf properties, which will help to discriminate 
between stellar and substellar models of the companion. 
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